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Abstract 

We present an extension of the window flow control analysis by R. Agrawal et. al. (Reference 
HI), C.-S. Chang (Reference jh]), and C.-S. Chang et. al. (Reference |[8l) to a system with random 
service time and fixed feedback delay. We consider two network service models. In the first model, 
the network service process itself has no time correlations. The second model addresses a two-state 
Markov-modulated service. 


I. Introduction 

The discovery that deterministic feedback systems can be expressed in the network calculus 
as the solution of a fixed-point equation in the dioid algebra was exploited for an analysis of 
a general feedback system [HI, O, [|8l. The analysis was conducted for deterministic systems, 
which cannot exploit statistical multiplexing of traffic flows. A detailed feedback model of the 
TCP Tahoe and TCP Reno algorithms as max-plus linear systems was presented in [|3l. The 
analysis accounts for randomness in the feedback system, but assumes that the underlying service 
process is deterministic. A further simplification is that the derivations in [[3]| are conducted for 
an overloaded system, with a claim that the results extend to any load condition. Overcoming the 
limitations of a deterministic analysis has motivated the development of the stochastic network 
calculus, where traffic and service are characterized by random processes [QAII . Whereas, over 
the last ten years, many key results of the network calculus have been shown to also hold 
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< w bits 



(delay d) 


Fig. 1: Network with window flow control. Traffic is admitted to the network only if its backlog 
does not exceed w bits. Feedback, with delay d, informs the network ingress about departures. 
Random processes A, A', D, and S, respectively, describe the external arrivals, the admitted 
arrivals, the departures, and the network service. 


under probabilistic assumptions, there has not been a successful treatment of stochastic feedback 
systems, see ifTTl p. 82], ||l2l p. 104], [[T3l pp. 215-216]. 

In this paper, we present an analysis of stochastic feedback systems with methods of the 
stochastic network calculus. Our paper extends the deterministic network calculus analysis of 
feedback systems in [[I|, [[6l, [[8l, llTSll . Following prior works, we study a window flow control 
system as shown in Fig. which enforces that only a limited number of w bits, referred to as 
window size, can be in transit at any time. The feedback consists of (acknowledgement) messages, 
which relay information on the amount of departing traffic back to the network entrance. The 
delay of the feedback is assumed to be d time units. The feedback signal opens or closes a 
throttle that prevents traffic from entering the network. 

Selecting a discrete-time model for the analysis, we derive upper and lower bounds on the 
service experienced by a traffic flow where the external arrivals and the available network service 
are characterized by bivariate random processes. We are interested in gaining insight into the 
service impediment caused by the feedback mechanism with given window size w and feedback 
delay d. This is done by deriving an equivalent service process, which describes the available 
service when taking into account the window flow control constraints. 

Our analysis applies the moment-generating function (MGF) network calculus from [|3- One 
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challenge for the analysis of a stoehastie feedbaek system with the MGF network ealeulus is 
that the standard method to eompute the frequently eneountered eonvolution operation is not 
well-suited for feedbaek systems. We adapt the MGF network ealeulus to feedbaek systems by 
deriving a sharpened version of the eonvolution estimate from [13 • 

It remains open to whieh degree our results ean be applied to other feedbaek systems. Even 
though our paper only presents single node results, we emphasize that available teehniques of the 
MGF network ealeulus permit an immediate extension to a multi-node network eonsisting of a 
sequenee of feedbaek systems. More eomplex network topologies, in partieular, nested feedbaek 
systems are not eovered. 

The rest of the paper is struetured as follows. In See. we provide baekground on network 
ealeulus, analysis of feedbaek systems, and the MGF network ealeulus with bivariate proeesses. 
In See.|^ we diseuss the obstaeles to a network ealeulus analysis of stoehastie feedbaek systems, 
and present the main results of this paper. In See. we eonsider a feedbaek system with an 
i.i.d. serviee model. See. [V| extends the analysis to a serviee model with time eorrelations. We 


present brief eonelusions in See. VI 


II. Bivariate Network Calculus 

The derivations in this paper will be done in the eontext of a bivariate network ealeulus, whieh 
has been used for deterministie as well as stoehastie analyses of networks [|71. We eonsider a 
diserete-time domain with f = 0,1, 2,... deseribing time slots. 

The vast majority of network ealeulus researeh is done with univariate funetions f{t) that 
quantify events in a time interval [0,f). A bivariate network ealeulus uses funetions f{s,t) to 
eharaeterize events in the time interval [s,t). The rationale for preferring the univariate network 
ealeulus is that it has stronger algebraie properties. On the other hand, univariate funetions do not 
lend themselves easily to probabilistie extensions. For example, eonsider a deterministie system 
that offers a time-invariant serviee S'(s, t) for eaeh time interval [s, t). Sinee S{s, t) = S{0,t — s), 
the serviee ean be eompletely deseribed by a univariate funetion S'(r) expressing the serviee in 
any time interval of length r. In the eorresponding probabilistie model, if the serviee is stationary, 
S{s,t) and S'(0,f — s) are equal in distribution, however, generally S{s,t) ^ 5(0, f — s), thus 
prohibiting the eonvenient reduetion to univariate funetions. 
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A. (A, ®)-dioid algebra for bivariate functions 

We denote by IF the family of bivariate functions that are non-negative and non-decreasing in 
the second argument. We use Fo to denote the set of functions f ^ F with /(f, f) = 0, which 
we refer to as causal functions. For f,g^F, the minimum (A) and convolution (®) operations 
are defined by 

f A g{s,t) = min{/(s, t),g{s, t)} , 
f (^g{s,t) = min {/(s, r) + g{T, t)} . 

S<T<t 

The two operations form a dioid algebra on the sets F and Fo [|7]|. The convolution operation in 
these dioids is associative, but not commutative. (In contrast, the convolution in the corresponding 
dioid algebra for univariate functions is commutative.) Also, the convolution distributes over the 
minimum, i.e., for three functions /, g, and h we have f ®{gf\h) = A/®/i{^The function 

{ 0 s > t, 

oo s < t , 

is the neutral element of the convolution operation, with f(S)6 = 6®f = f for any f E F. 

For the analysis of feedback systems, we need to convolve a function multiple times with 
itself. For that, we use the notation 

= /'” = /, 

f{n+l) ^ j{n) ^ y ^ (n > 1) . 

We define the subadditive closure of f E F, denoted by /*, as 

OO 

/* = 5a/a/(2)a/(3)a... = /\/("). 

n=0 

The attribute ‘subadditive’ derives from the property that f*{s,t) < /*('S,r) -f f*{T,t) for 
s < T < t. Every subadditive function / satisfies / <8)/ = /.If/€ Fo, then its subadditive 
closure is given by f* = lim n^oo f^ 

'statements that list functions without arguments hold for all pairs (s, f) with s < t. Also, we give the convolution operation 
precedence over the minimum, which allows us to omit some parentheses. 
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B. Bivariate arrival and service processes 

The arrivals to a network element are eharaeterized by a bivariate proeess A, where 74 (s, t) 
deseribes the eumulative arrivals in the time interval [s,t). We use D(s,t) to deseribe the 
departures in the time interval [s,t), subject to the causality condition D{0,t) < Both 

A and D are causal functions {A, D G Fo)- If Ofc and dk denote the arrivals and departures, 
respectively, at time k, we have 

t-i t-i 

A{s,t) = ''^cik and D{s, t) = dk . 

k=s k=s 

The available service at a network element is described by a bivariate service process S'(s, t) G 
Jo that satisfies the input-output relationship 

D>A(^S 

for any arrival process and corresponding departure process at that element. Such a process 
is called a dynamic server in 0 p. 178]. An exact service process satisfies the input-output 
relationship with equality. The service offered by a sequence of network elements with service 
processes Si, S 2 , ■ ■ ■, Sn is given by their convolution S'! ® S '2 ® ® Sn. 

Given arrival and service processes, one can formulate bounds on backlog at a network 
element. The backlog 5 of a network element denotes the arrivals that have not yet departed, 
given by B(t) = A(0,f) — D{0,t). A bound on the backlog can be expressed in terms of the 
deconvolution operation (0), which, for two functions /, (7 G J”, is defined by 

/ 0 g{s, t) = max {/(r, t) - g{T, s)} . 

0<T<S 

With this operation we have [13 Theorem 2] 

Bit) <A(Z)S(t,t). (1) 

Bounds on the delay and the burstiness of departures also involve the deconvolution operation. 

C. Bivariate feedback systems 

A network element with feedback is one where the departures from the element influence the 
arrivals to the same element. Fig. depicts a generic system-theoretic model of a closed-loop 
feedback system. There is a network element with service process F whose output is re-combined 
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Fig. 2: Generic feedback server. 


with the external input, so that the arrivals to the network element are the minimum of its output 
and the external arrivals. Such a system will be referred to as feedback server. The feedback 
server consists of a (non-causal) service element with service process F{s,t) (F ^ Fo) and 
arrival and departure functions D' > D ® F. The service element labeled by ‘min’ in Fig. 
represents a throttle, which enforces the minimum D = AA D'. With this, A, D, and F satisfy 
D>AAD(^F. As shown in [O, [|3 for univariate functions, the closed-loop system can be 
replaced by an equivalent system, which consists of a single network element without feedback 
with service process F*, where F* is the subadditive closure of F. In dHl, the result has been 
extended to bivariate functions, as expressed in the following lemma. 

Lemma 1: (see /l^ Lemma 2.2]) Given a feedback server as in Fig. with bivariate functions 
A, D, and F. If F ^ Fo and D > A A D ® F, then D > A® F*. 

Since the lemma applies when A, D, and F are random processes, an analysis of stochastic 
feedback system appears readily available. However, computing the subadditive closure F* for 
a bivariate random service process presents considerable difficulties. 

D. MGF network calculus 

The MGF network calculus ||9l offers an analysis of network elements, when the arrivals and 
the offered service are bivariate random processes that are characterized in terms of their moment¬ 
generating functions. The MGF network calculus has been frequently applied to the analysis of 
wireless networks, e.g., |l2l, ifTO . lITbll . [UHl, since the random service model can capture the 
randomness of a wireless transmission system. We denote the moment-generating function of 
a random variable X for any 6* G M by Mx{d) = E[e^^]. The MGF calculus exploits that for 
any two independent random variables X and Y, the relation Mx+y{0) = Mx{0)My{0) holds. 
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The stochastic analysis of feedback systems in this paper will take an MGF network calculus 
approach. 

The moment-generating functions of an arrival process A and a service process S' for 6 * > 0 
are denoted by 

Ma(0, s, t) = and Ms{-e, s, t) = . 

We assume that the arrival and the service are independent. For characterizing the random service 
of a feedback system, we will use a bivariate version of the statistical service curve from [[5l 
and the effective capacity from [|T71 . A statistical service curve of a bivariate service process 
S for a given e > 0 dH is defined by the property that 

Pr ^S(s, f) < 5^(s, f) j < £ . (2) 

The statistical service curve is a deterministic function giving a lower bound on the available 
service that is violated with a probability e or less. Using the Chemoff bound, a statistical service 
curve can be computed from Ms as 

S%s,t) = max^jloge - log I . (3) 

An alternative measure to describe a stochastic service in terms of Ms is the effective capacity 
'ys{—d) IHtII . defined for 6 * > 0 by 

Isi-O) = -^\ogMs{-0,OA) ■ (4) 

t^oo 

Since it is defined as a time limit, the effective capacity is most useful when reasoning about 
long-term traffic rates and scaling properties. In particular, 75 ( 0 ) = lim 5 i_^o 7s(“^) equals the 
average service rate. In contrast, the statistical service curve provides a bound for finite values 
of t. 

Given the moment-generating functions of two bivariate processes / and g, bounds on the 
convolution and deconvolution, given in Il9l. are 

t 

Mf^g{-0, S,t) Mf{-e, s, T)Mg{-e, T, t ), (5) 

T=S 

s 

Mf0g{9, s,t) Mf{9, T, t)Mg{-9, T, s). (6) 

r=0 
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The convolution bound is used to estimate the moment-generating function of a sequence of 
service elements. The deconvolution bound plays a role when computing performance bounds. 
For example, with an application of the Chernoff bound, we can obtain from Eq. o a bound 
on the backlog distribution Pr{B{t) > b*(t)) < e [|9||, where 

jlog - logej . (7) 

The delay can be treated in a similar fashion. 

III. Towards a Stochastic Feedback Analysis 

The generalization of an analysis of a deterministic feedback system to random arrival and 
service processes has remained open for considerable time. In this section, we present results 
that make such an analysis possible. We describe the issues that make a stochastic analysis of 
feedback systems within the framework of the network calculus hard, and then address how to 
resolve them. 


A. Model Description 

We will analyze the network with window flow control in Fig. Traffic with arrival process 
A{s,t) is serviced by a network element with service process S{s,t), subject to the additional 
constraint that the total backlog in the element at any time may not exceed w > 0. Traffic in 
excess of that constraint is held in a FIFO buffer at the network entrance. The traffic leaving the 
network is expressed by a departure function D(s,t). The feedback information consists of the 
value of the departure function delayed by d > 0 time slots, i.e., D{s,t — d). If we use A' to 
denote the arrived traffic that is admitted into the network, the flow control system requires that 

A'(s, t) = min {A(s, f), D{s, t — d) + w} . (8) 


This control ensures that admitted traffic that has not yet departed cannot exceed w. 


The system can be described by a feedback model as discussed in Subsec. II-C Let the 
network service be given by a service process S, i.e.. 


D> A'® S. 


(9) 
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win 



Fig. 3: Model with flow control with window size te > 0 and feedback delay d > 0. 


Further, for te > 0 we define a function <5+"' by 




{ w s > t, 

oo s < t . 


( 10 ) 


This allows us to write f{s,t) + te = / ® f). Note that 5+"' is not a causal function, 

i.e., > 0 for t < s. Even though the convolution of bivariate functions is generally not 

commutative, we have / ® = 5’'""' (8) / for a bivariate process /. Lastly, a service element 

offering a delay of d > 0 is given by a service process 5d = S{s,t — d), such that f{s,t — d) = 
f 8 dd{s,t) for every / G With these definitions we rewrite Eq. as 


A = AAD®5d®5+'^. 


( 11 ) 


In Fig. 1^ we illustrate the feedback system as a concatenation of service elements. By inserting 
Eq. @ into Eq. ( [TT] ) we obtain 

A'> A^A'® S . (12) 

Note that A' satisfies the conditions of Eemma[^ with F = S ® dd® d^'^. Applying the lemma 
and inserting the result into Eq. @ yields 

D>A®{S®dd® d+^y ® S . 


Therefore, the flow control system can be represented by an equivalent service process S'„in 
given by 

S^in = {s ®dd®d^^y ® S . (13) 


If S is an exact service process, so is S^n 
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B. Challenges of stochastic feedback systems 

The analysis of a feedback system with a bivariate random service process S requires the 
computation of the subadditive closure in Eq. ( [T3] ). Let us consider for the moment a system with 
no feedback delay, that is, d = 0. Since {S ® 5 +"') ® (S' ® < 5 = S' ® S' ® the computation 
of the subadditive closure involves the convolution of S with itself. When the service process 
is subadditive, we have S' = S' ® S, and therefore obtain (S' ® ^ S = S. Hence, when 

a subadditive process S' is the service process of a closed-loop system without feedback delay, 
feedback has no impact on the overall service. On the other hand, a service process with a 
nonzero feedback delay is generally not subadditive and, therefore, the convolution of such a 
process will not yield a trivial result. To see that the concatenation of S' and 5d with d > 0 is 
not subadditive, we compute 

S ® Sd{s, t) + S ® 6d{T, t) = S{s,T - d) + S'(r, t - d) 

< S{s,t — d) 

= S' ® Sd{s,t) . 


To avoid trivial cases, we will henceforth consider feedback systems with d > 0. 

Writing the expression for the equivalent service process in Eq. ( [T3] ) as 

OO 

S'win = /\ 

n=0 

makes apparent the need for tools to analyze the distribution of minima involving a bivariate 
random service process S'. Note that the n-fold convolutions themselves involve minima as well 
as sums of random processes. 

Consider, for comparison, the same formula in the case of an univariate service process S'(f). 
There, one can exploit the commutativity of the convolution to obtain for the n-th term in the 
minimum 



(^ ® dd ® ® S{t) = 0 0 6nd{t) = - nd) + nw . 

In particular, if S' itself is subadditive, then S'^”^ = S', the n-th term is S{t — nd) + nw, and the 
entire minimum reduces to S^in = S ® So, where So{t) = w\^ - Since deterministic feedback 
systems can generally be expressed using univariate service processes, the computation of the 
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subadditive closure of the service process is much simplified. However, for bivariate service 
processes, the convolution is not commutative, and we must compute the n-th term as 


(S'(g) 0 0 S'(s,f) = min ^S(ri_i, r* - d) j + S(rn, f) | + nte , (15) 

where the minimum ranges over all non-decreasing sequences To,... ,Tn with To = s and Tn <t. 
One issue with this expression is that the number of terms grows rapidly with n and t — s. 
The second difficulty is that the joint distribution of S'(ri_i,rj — d) for i = 1,... ,n is not 
determined by the distribution of S{s,t) alone, but requires information on time correlations. 
We now proceed to address these problems. 


C. Service processes of stochastic feedback systems 

Our first result provides an exact characterization of S^m for the special case d = 1. 


Lemma 2: Consider a feedback system as in Fig. with an additive service process 

t-i 

S{s,t) = ^Ck, 

k=s 

where {ck)k>i is an arbitrary sequence of non-negative random variables. For d = 1 and te > 0, 
the equivalent service process is given by 

t-i 

5'win(s, f) = ^ niin{cfc, te} . 

k=s 


Proof: Set d = 1. For n = 1, the convolution in Eq. ( [15] ) is given by 

S' g) (5i g) 5^'^ ® S{s, t) = min {S{s, r — 1) -f S(r, f)} -f w 

S<T<t^ ^ 

(T-2 t-l s 

“ k = S ^ 


= mm 

S<T<t 


+ w 


k=T 


= S(s, t) — max Ck + w 

k£[s,t) 


Here, the convolution replaces the largest value in the sum representing S{s,t) by w. Likewise, 
for n > 1, the minimum in Eq. ( [T5| ) is obtained by replacing the n largest values of Ck on [s, t) 
with w. In this case, the minimum in Eq. is attained for n = ^{k E [s,f) | Ck > w}. ■ 
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The following theorem shows that the exact result from Lemma provides a lower bound on 
Swin for d > 1 in terms of the ratio w/d. The theorem also provides a complementary upper 
bound on S'win- 

Theorem 1: Given a service process S G J^o- Let be the equivalent service process for a 
feedback system with parameters d > 0 and te > 0. Then 

S'{s,t) < S'win(s,f) < min 1"^] te} , 0<s<f, (16) 

where S' is the equivalent service process (given by Eq. ( [T?] )) with d' = 1 and w' = w/d. 

Proof: For the lower bound, observe that 

= (d, (17) 

This implies 

r 

S'(g) drf 0 d+"'(s,f) > min < fs(rj_i, r* - 1) + w') 

S=To<"-<T^=t I V / 

^ i=l 

/ l\W 

= (s,f). 

Indeed, the left-hand side appears as a term in the minimum on the right-hand side. Therefore, 
the n-th term in the expression for S'win in Eq. ( [T?] ) is bounded from below by the (nd)-th term 
in the expression for S'. 

For the upper bound, we simply use just two terms from the minimum in Eq. ( fT4| ) to bound 
S'win, namely n = 0 and n = \. ■ 

An implication of the theorem is that if S', S'win, and S'Lj^ have long-term average rates C, 
Gwin, and GL„, then 

G^in < (S'win < minjG, . 

Theorem holds for general (deterministic or random) service processes. The bounds on the 
average rates depend only on the ratio w/d, not on the values of the individual parameters. 
Our main results, in Theorems and will strengthen the lower bound when d > 1 for two 
important classes of additive service processes. They rely on the next lemma, which reduces the 
number of terms that contribute to Eq. ( [T?] ) and limits the range of the minimum in Eq. ( [T4| ) to 

"< M- 
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Fig. 4: Geometric interpretation of the index set Cn{s,t). 


Lemma 3: Given a feedback system with a service process S G 3Fo- Then, for any choice of 
d > 0 and te > 0, 


5'win(s,^)= /\ 


n=0 


min S{Ti_i,Ti -d) + S{Tn,t) + 

\ *=1 


nw 


(18) 


where the minimum in the braces ranges over 



T'oi 'Tl, 


I 'Tn. 


Tn 


s, Tn < t, and Ti — Ti-i > d for i 





if nd < t — s. \i nd > t — s, then Cn{s,t) contains the single sequence r* = (s + id) A t for 
i = 0,... ,n. 


The set G„(s, t) has a geometric interpretation, illustrated in Fig. Each sequence To,... ,Tn 
in Cn{s,t) corresponds to a collection of n disjoint subintervals of length d in [s,t], given by 
(tj — d, Ti]. On each of these subintervals, the original service process S'(rj — d, Ti] is interrupted 
by a delay of length d, followed by an addition of w. In the special case d = 1, the set G„(s, t) 
consists precisely of the n-element subsets of [s,t), and we recover the statement of Lemma 

Proof: We need to prove that only sequences Tq, ..., in G„(s, f) contribute to the minimum 
in Eq. \f nd > t — s, the minimum value is achieved by subdividing [s, f) into n subintervals 
of length at most d, so that all terms involving S vanish. Then the right-hand side of Eq. 0 
equals nw > \w. 

Eor nd < t — s, we proceed by induction over n. The n = 0 term in both Eq. ( [T?] ) and Eq. ( [TS] ) 
is given by S{s,t) and there is nothing to show. The n = 1 term equals 

S ^ 6d® ® S{s, t) = min {S{s,t — d) -f S'(r, f)| + w . (19) 

The key observation is that the expression in the braces is non-increasing in r for s < t < d, 
since the first term vanishes while the second term is always non-increasing. Therefore, the 
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minimum is achieved for some t > s + d. This means that the sequence Tq = s, Ti = r lies in 


Now consider n> 1. Suppose we already know that 


(S (g) 6d® ® S(s,t) = min 

To,...,T„&Cn(s 


n 

j -d) + + 

i=0 ^ 


nw 


for all t > s. Using the associativity of the convolution, we write 


{S®5d® ®S=\{S^5d^ 


and then expand 


{S(g6d<g6+^Y^+^^ (gS{s,t) 


min 

s<r<t 


iSg)6dg> ® (s, T-d) + S{t, 


+ w. 


For T — d< s + nd, the term in the square brackets takes the constant value nw, while the second 
summand is non-increasing in r. Therefore, the minimum occurs at some r > s + {n+ l)d. By 
the inductive assumption. 




d) 


min 

To,---,Tn&Cn{s,T — d) 


- d) + S{Tn,T 

i=\ 



Since Tq, ... 


rn,r 


e Cn+i{s,t), the claim is proved. 


IV. Variable Bit Rate Service with Feedback 

We next apply the results from the previous section to a specific random service process in 
a feedback system, consisting of a work-conserving FIFO buffer with a random time-variable 
service rate, which we refer to as variable bit rate (VBR) server. The feedback mechanism is 
as described earlier with window size tc > 0 and feedback delay d > 0. The VBR server offers 
the service process S{s,t) = where Ck is a random amount of available service in 

the k-\h time slot. We assume that the c^’s are independent and identically distributed (i.i.d.) 
random variables, with moment-generating functions Mc{9) = The moment-generating 

function of S is Ms{0,s,t) = {Mc{0)Y~^, and the effective capacity has the simple expression 
gs{-0) = -|logMc(-d). 
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A. Bounds for a VBR service with feedback 

We now derive bounds on the equivalent serviee proeess Swin of a VBR server in the feedbaek 
system of Fig. The bounds will be expressed either in terms of a statistieal serviee curve for 
a chosen violation probability e > 0 (using Eq. @) or in terms of the effective capacity (from 
Eq. Q), both of which require bounds on the moment-generating function of S^i^, denoted by 
^Swin- We use the expression for given in Eq. ( [T4| ) as the starting point for the computation 
of where we exploit the following relationship of moment-generating functions. 

Lemma 4: Given two bivariate random processes / and g. Then, for every 6^ > 0, 

MfAgi-0, s, t) < s, t) + Mg{-9, s, f). 

Note that this lemma does not require / and g to be independent. 

Proof: 

= Mf{-9, s, t) + Mg{-9, s, t). 


Since the service guaranteed by S is independent on disjoint intervals, that is S(ti,t 2 ) and 
S{t 3 ,ti) are independent for h < t 2 < < ti, standard techniques for bounding the moment¬ 

generating functions of sums and convolutions can be applied directly to Eq. ( fT?] ) to obtain, for 
9>0, 


Me 


d)y Ms{-9,Tr„t)t 

n = 0 S=To'^...<Tn'^t \ 2=1 

< (M,(-9))‘- f; (‘ ■ " t ^ ^ ((W(-9))-V''")’ 

n=0 

(M,(-0))‘-* 


—Onw 


n 


( 20 ) 


(1 - (Mc(-0))-'^e-®“’)‘“'*+^ ’ 

as long as the convergence condition 


< 1 
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holds. The first line follows by Lemma|^ and the independence of the service on disjoint intervals 
in each summand. In the second line, Ms is expressed in terms of Me, the binomial coefficient 
counts the number of non-decreasing sequences Tq < ... < Tn with To = s and Tn <t. The last 
line follows from the identity 


OO 

E 

n=0 


t + n 
n 


X = 


(1 


X 


it+l 


(|x| < 1). 


Eq. pOl ) provides a useful bound only when 9 is chosen so that the convergence condition is 
satisfied, and the right-hand side of Eq. pO] ) is less than one. We now improve the bound with 
the help of Eemma 


Theorem 2: Eet S{s,t) be a VBR server with feedback as described above. Then, for every 

e>o, 

Ms^J-e,s,t) < + . ( 21 ) 


Proof: Eix 9 > 0. Since S„^ia{s,t) is non-decreasing in t, we can round down the length of 
the time interval to the nearest integer multiple of d. Thus, it suffices to consider the case where 
t — s = Nd for some integer N. By Eemma 


N 


5'win(s,f) = /\ 


n=0 


mm 

To,...,rn£Cn{s,t) 


n 

-d) + S{Tn,t)^ +nw 

i=l 


By Eemma and using that S is independent on disjoint time intervals, we estimate 


N 


Ms,,;^{-9,s,t) <^< ^ iY[Ms{-9,Ti-i,Ti-d)]Ms{-9,Tn,t)t 


—Onw 


^=0 ro,...,rneCn(s,t) 

The number of sequences in Cn{s,t) is bounded by 


\Cn{s,t)\ = 


(N — n)d + n 


n 


j-j- {N — n)d -f j ^ j-j- (iV — n)d + dj _ ^^ 


i=i 




Expressing Ms in terms of M^, we arrive at 


N 


n=0 

n=0 ^ ' 

= [{Me{-9)f ^ de-^^Y . 


16 






Since N = the elaim is proved. 


Theorem 1^ (as well as Eq. (|20|)) ean be direetly inserted into expressions of the MGF ealeulus, 
e.g., for the baeklog expression in Eq. Q. Further, via Eq. Q, the theorem also provides a 
statistieal serviee eurve for S'„in, whieh we will denote as By taking the logarithm of 
M^win, we ean obtain bounds on the effeetive eapaeity 7 win, as expressed in this eorollary. 

Corollary 1: Eower bounds on the effeetive eapaeity 7 ^^ of a VBR proeess S{s,t) with 
feedbaek are given for 9 > Oby 

a) 7win(-0) > 7s(-^) + I log (1 - , (22) 

b) 7win(-0) > 7 s(-^) - ^ log (1 + . (23) 


The first bound follows from Eq. ( |20| ), and requires the eonvergenee eondition 'js{—9) < 
The seeond bound follows immediately from Theorem]^ Eqs. and ( |2^ elearly express the 
serviee impediment due to the feedbaek proeess, by subtraeting a positive term from the available 


serviee without feedbaek. In Subsee. IV-C we present numerieal examples that evaluate both 
bounds. 

All results in this seetion ean be applied to a ‘leftover serviee’ model at a server offering a 
eonstant-rate serviee in the presenee of eross-traffie with independent inerements. The leftover 
serviee expresses the serviee available to a flow in terms of the eapaeity that is left unused 
by eompeting eross-traffie. The leftover model assumes that the analyzed traffle flow has lower 
priority than the eross-traffie. Explieitly, if the serviee rate is C and eross-traffie arrivals are 
given by 

t-i 

A%s,t) = Y.' 




k=s 


where the eross-traffie arrivals in eaeh time slot are given by an i.i.d. sequenee of random 
variables a^, then the leftover serviee available to the flow satisfies 

t-i 

>Y,(C - al) . 

k=s 

Although the summands C — al may take negative values. Theorems and remain valid and 
provide non-trivial bounds on the serviee proeess with feedbaek, as long as the stability eondition 
E[al] < G is satisfied. 
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B. Quality of the VBR bounds 

We next use Theorem[^to address the aeeuraey of the derived lower bounds from Subsec. IV-A 
Applying the theorem to the VBR server gives the bounds 




t-l 


k=s 


W 


^minjcfc, ^ \ < S'win(s,f) < min 


in| Vcfc, 


k=s 


d 


w 


(24) 


Even though the upper bound is optimistic, its difference to the (lower) bounds on the service 
computed from Theorem limits the deviation from the true value of S'win- Moreover, the 


difference between the upper and lower bounds in Eq. (24) indicates the range of useful estimates. 


Note that the lower bound in Eq. ( |24l ) is exact for the special case d = 1. In particular, any 
lower bound that falls below this bound can be replaced by the simpler estimate of Eq. ( |24| ). 

Eet us, for the moment, consider a deterministic server with Ck = C. Then, the lower and 
upper bounds are essentially equivalent, giving S'„in(s,f) min{(7, — s). This corresponds 

to a well-known expression for the throughput of a window flow control system, e.g., dH Eq. 
6.1]. On the other hand, when S is random, the upper and lower bounds have different long-term 
average rates. 

We will use Eq. ( |24| ) to obtain upper and lower bounds on and 7„in. We obtain a lower 
bound on by taking the moment-generating function of the lower bound in Eq. ( |24| ), resulting 
for 6* > 0 in 

(25) 

Note the change of direction of the inequality since we use —9 as a function argument. Inserting 
this bound into Eq. ([^ provides a lower bound on 5^;j^(s,f). 

Obtaining an opposing upper bound on from Eq. ( |24l ) is not as straightforward, since 
statistical service curves express lower bounds on the available service. We exploit that, for 
some VBR servers, it is possible to get the exact distribution of S{s,t) = where 

S corresponds to the service process of the VBR system without feedback. Eor example, in 


Subsec. IV-C where we use an exponentially distributed Ck, the service S{s,t) has an Erlang 
distribution. In such cases, since the upper bound in Eq. ( |24| ) implies S'win < S', the e -quantile 
of S is an upper bound for any statistical service curve of the VBR server with feedback. 


Thus, Eq. (24) gives 


‘5win(s,t) < min I £-quantile of S'(s,f), 


d 


w 




(26) 
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where we also use the seeond (deterministie) term of the upper bound in Eq. ( [24| ). 

Bounds for the effeetive eapaeity are direetly obtained from Eq. ( [24| ) by first eomputing 
moment-generating funetions and then taking the logarithm. 

Corollary 2: Under the assumptions of Theorem]^ 7 win(—6^) is bounded for 9 > Oby 

<7wl„(-9)<min{7s(-9),^} . 


The lower bound beeomes an equality when d = 1. The upper bound is strietly larger than the 
lower bound for all values of 0 > 0. Most importantly, the upper bound beeomes sharp if we 
take d —)■ cx) while holding w/d fixed. In the deterministie ease, with Ck = C, the upper and 
lower bounds both equal min{U, ^}, in aeeordanee with our previous diseussion. Moreover, in 
the limit 6 ^ oo, the right-hand side of Eq. ( |2^ is asymptotie to the shown bounds. Therefore, 
the bounds on the effeetive eapaeity are sharp in the limit 6* —)■ cx). 


C. Numerical evaluation of VBR bounds 

We now present a numerieal evaluation of our bounds for the VBR server with feedbaek. We 
eonsider a VBR server with an exponential distribution with moment-generating funetion and 
effeetive eapaeity given by 

Msi-e, s, f) = (1 + , ysi-0) = ^ iog(i + ce). 

We seleet C = 1 Mb for the available serviee in a time slot of length 1 ms, whieh gives an 
average serviee rate of 1 Gbps. 


In Eigs. 5(a) and 5(b) we plot statistieal serviee eurves as funetions of time, where 

we use a violation probability of e = 10“®. We eompute multiple serviee eurves where we vary 
the delay d and the window size w, but fix the ratio w/d. The statistieal serviee eurves use 
Eq. Q with the bound on from Theorem We eompare these serviee eurves to the 

upper and lower bounds eomputed from Eq. 


Eig. 5(a) shows statistieal serviee eurves, plotted as solid lines, using a ratio w/d = 100 Mbps. 
The upper and lower bounds are represented by dash-dotted lines. We also indieate the two terms 
in the upper bound of Eq. ( |2^ by dotted lines. Note that the lower bound is in faet equal to 
for the speeial ease d = 1 ms, due to Eemma 

It is evident that the rates of the serviee eurves mateh well to those of the upper and lower 
bounds. We observe that inereasing the delay and the window size simultaneously improves the 
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Time (ms) 


(a) w/d = 100 Mbps. 



Fig. 5: Statistical service eurves for VBR server with feedbaek (Avg. rate: 1 Gbps, 

e = 10~®). The dotted lines correspond to the two terms of the upper bound in Eq. (|2^. 


available service. For large values of d, the service eurves appear to have a supremum well 
below the plotted upper bound. This is expected, sinee the rates of our bounds are exact for 
d —)■ oo. In a deterministic feedback system with a fixed rate of 1 Gbps, varying w and d with 
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(a) w/d = 100 Mbps. 



Fig. 6: Effective capacity 7win(6*) of VBR server with feedback (Avg. rate: 1 Gbps). 


a fixed ratio w/d results in all cases in an essentially constant rate service of 100 Mbps, with 
minuscule deviations. 


Fig. |5(b) evaluates the same scenario for a different parameter selection, this time, fixing 
w/d = 500 Mbps. Note that the range of the y-axis is modified from Fig. 5(a) For the shown 
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range of time values, the rate of the statistieal serviee eurves are elose to that of the lower bound, 
but notieeably smaller than that of the upper bound. This will ehange when we eonsider larger 
time intervals. The reason is that with this ehoiee of w/d, the upper bound of the serviee eurves 


is dominated for a longer period of time by the first term of Eq. ( [26| ). Onee the seeond term 
(with rate w/d) governs the bound, the rate of the statistieal serviee eurves will be elose to w/d 
as well. 

We now turn to the effeetive eapaeity. For the same set of parameters as before, we evaluate in 
Fig. I^the lower bounds of the effeetive eapaeity 7win from Corollary Note that the eorollary 
presents two bounds. Corollary a) (Eq. ([22l)) does not take advantage of Theorem whereas 
Corollary [^b) (Eq. ( |2^ ) involves Theorem]^ We plot these bounds on 7win(—as a funetion of 
0 > 0. Reeall that the aetual (not estimated) effeetive eapaeity 7win(—6*) is a deereasing funetion 
of 9, and that its value for 0 —)■ 0 is the average serviee rate. 


Fig. 6(a) depiets the bounds for w/d = 100 Mbps, and Fig. 6(b) those for w/d = 500 Mbps. 
Consider first the lower and upper bounds from Corollary whieh are indieated by dash-dotted 


lines. As diseussed in See. IV-B, the bounds eonverge for large values of 9. The figure indieates 
that the eonvergenee oeeurs early. We observe that bounds from Corollary [^a) are inferior to 
those of Corollary [^b), whieh emphasizes the value of applying Theorem]^ For d > 1 ms, the 
lower bounds of Corollary [TJb) for different values of w and d are aeeurate when 9 > A - 10“^ 
(Fig. 6(a) ) and 6* > 0.5 ■ 10“® (Fig. |6(b) ), but degrade for small values of 9. When this happens 
the lower bound from Corollary should be used. For d = 1 ms, the bounds from Corollary [TJb) 
are pessimistie over a large range. Here, the lower bound from Corollary whieh is exaet when 
d = 1 ms, is the better result. Sinee the aetual 7win(—inereases when redueing 9, the best 
estimates of the average available serviee of are obtained at the maximum of the eurves. 
It is interesting to observe that Corollary generates bounds whieh are not always monotonie 
when inereasing d. 

Our last numerieal example presents a baeklog analysis of the VBR server with feedbaek, 
by applying Eq. 0. For the arrivals, we seleet a proeess similar to the serviee proeess, where 
arrivals in eaeh time slot are i.i.d. with an exponential distribution that generates arrivals at an 
average rate of A Mbps. The moment-generating funetion of this arrival proeess, denoted by 
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Fig. 7: Backlog bounds for VBR server with feedbaek {C = 1 Gbps). 


Ma, is given by 

The serviee is the same VBR server as evaluated before. By setting <7 = 1 ms, we ean take 
advantage of the exaet expression for that follows from Lemma Even though we have 

exaet moment-generating funetions for arrivals and serviee, the baeklog bound of the MGF 
network ealeulus will be pessimistie, due to the deconvolution expression in Eq. In Eig. |7] 
we depiet baeklog bounds for the ratios w/d = 100 Mbps and 500 Mbps, for different violation 
probabilities e, as funetions of the arrival rate A. We note that this presents the first probabilistie 
baeklog bounds of a random feedbaek system using methods of the network ealeulus. We observe 
that for eaeh ehoiee of w, the system saturates at a well-defined rate. The saturation rates are 
elose to the lower bounds on 7win(—6^) for 0 —)■ 0 in Fig. The plots indieate a low sensitivity 
of the bounds to £. We also inelude simulation results for £ = A comparison with the 

simulations shows that the aehieved baeklog bounds are pessimistie, but traek the blow up of the 
baeklog at high utilizations well. The analytieal bounds will be more pessimistie for e = 10“^ 
and less so for e = 10“®. 

^The simulations represents runs of 10® time slots, where the simulation is started with an empty system, and results from 
first 10® time slots are discarded. 
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V. Service Processes with Positive Time Correlations 


The VBR server from the previous seetion offered an i.i.d. serviee in eaeh time slot, and the 
time eorrelations of the feedbaek system resulted exelusively from the feedbaek meehanisms. 
Nonetheless, the analysis of this simple feedback system proved to be substantial. This raises the 
question whether feedback systems with more complex service processes are at all tractable with 
our analysis approach. In this section, we provide a positive answer, by analyzing a feedback 
system for a service process with memory. 

The server is represented by a Markov-modulated On-Off (MMOO) process [jTl, which oper¬ 
ates in two states. In the ON state (state I) the server is transmitting a constant amount of P > 0 
units of traffic per time slot. In the OFF state (state 0), the server does not transmit any traffic. 
The state is selected at the beginning of each time slot using fixed transition probabilities, where 
Pij denotes the probability of moving from state i to state j (i,j G {0,1}). By definition, the 
service process is additive, 

t-i 

S{s,t) = ^Cfc, 

k=s 

where = P if the system is in the ON state at time k, and = 0 otherwise. Since poi = 1 ~Poo 
and piQ = 1 —pii, the system is fully characterized by the three parameters, poo, pn, and P. The 
Markov chain is assumed to be in its steady-state, where the probability that the system is ON 
in any given time slot equals Pr{ON) = • MMOO processes are frequently used in the 

literature for modelling bursty traffic or service [fT4ll . The average rate of the MMOO process is 

E[S{s,t)] ^ Poi p 
t- s Poi+ Pio 

Its effective capacity is given in (71 Eq. ( 7 . 18 )] as 
7si-0) = Q {(Poo + ^/ (poo + Piie-®^)2 - 4(poo + Pii - ( 27 ) 

for all 0 > 0. 

We assume throughout this section that the transition probabilities satisfy poi +Pio < 1- This 
condition ensures that the system does not alternative rapidly between the two states. Smaller 
values of poi cause the system to linger in the OFF state, while smaller values of pio cause 
extended bursts. 


24 






A. Time-correlation properties of MMOO processes 

The analysis of an MMOO server with feedbaek requires a deeper inspeetion of the properties 
of an MMOO proeess. The following lemma eoncems the time correlations of the underlying 
two-state Markov chain. In this subsection only, we allow the parameter 9 to take both positive 
and negative values. 

Lemma 5: Consider a two-state Markov chain with poi +Pio < 1- For every strictly increasing 
sequence ti, ... ,Tn, the probability that the system is in the ON state at times ri,...,is a 
decreasing function of the differences r* — Tj-i, i = 2,... ,n. 


Proof: Consider the state variables 

{ 1 if the state is ON at time k, 

0 else. 

The product X^-. is the indicator of the event that the system is in the ON state at each time 
Ti, for i = 1,... ,n. By the Markov property, 


Pr (ON at time r* for all i 


l,...,n) = E 


n 



, 2=1 


= Pr(A'„ = 1) ■ n® = 1] ■ 

i=2 

By stationarity, the leading factor is a constant determined by the steady state of the Markov 
chain, and the i-th factor in the product depends on — rj_i. We now verify that each of these 
factors decreases with r* — n-i. 

Let p = Pr{Xk = 1) be the probability that the state is ON at time k, and let r = Ti — Ti-i. 
We compute the r-step transition matrix as 


Poo Poi 
Pio Pn 



p —p 

fl-p) (1-p) 


(28) 


where the first matrix on the right-hand side is the spectral projection onto the steady state, 
/i = 1 — poi — Pio is the non-trivial eigenvalue of the transition matrix, and the second matrix is 
the spectral projection onto the eigenstate corresponding to p. Since 0 < p < 1, 


E [X^^ I Xr^_^ = 1] = p + (1 - p)p^ 
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is a decreasing function of r, proving the claim. ■ 

The next lemma provides bounds on the moment-generating function of the MMOO service 
process. 


Lemma 6: Let S be an MMOO process as described above. If Poi+Pio < then the following 
inequalities hold for all 6^ G M: 

1) For every strictly increasing sequence ti < ■ ■ ■ < Tn, 


E 


< Ms(6»,0,n). 


2) Ms(6*,0,t) is supermultiplicative in t, 

Ms{9, 0, s) ■ Ms{9, 0, t) < Ms{9, 0,s + t) (Vs, f > 0) 

3) The moment-generating function is bounded by 

{M,{9)y < Ms{9,0,t) < (Vt > 0), 


(29) 


(30) 


(31) 


where Mc{9) = (1—p) +pe^^ is the moment-generating function of the service in a single 
time slot, and where m+(6*) is the larger eigenvalue of the matrix 


m = 


Poo Poi 
PlO Pll 



(32) 


4) Furthermore, 

Ms{9,0,t)>K{9){m+{9)Y, 

where 0 < iV(6*) < 1 is an explicit constant (to be computed in the proof). 


(33) 


Proof: 

1) Consider first the case 6* > 0, and let be the indicator function that the system is in 
the ON state at time k. Writing 

Ck = PXk, = l + {e^^ -l)Xk, 

we expand 

n 

o,, _ ^ ^ _ 1)1 J| ^ 

i=l j'eJ 
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For each subset J C {ri,... ,r„}, the distance between consecutive elements increases 
with the distances — rj_i for i = 2,... ,n, and all coefficients are positive. Therefore, 
we can apply Lemma to see that 


E 


M 

1 





76 J 


is a decreasing function of r* — rj_i. Since these differences take the smallest possible 
value when {ri,..., r^} = {0,..., n — 1}, this proves the claim for 6 > 0. 

For 6* < 0, let Yfc = 1 — be the indicator function that the state is OFF at time k, and 
set (j) = —6 > 0. We write 

expand the product as a sum 

ieJ 

and argue as in the other case. 

2) Fix 0 G M. For integers i > 0, let f{£) = E _ Factoring the 

exponential as in the proof of Eq. p9| ), it follows from Lemma that f{i) decreases 
with i. Taking £ —)■ oo and i = 0 yields Eq. ( [^ . 

3) By Part 2, the function g{t) = \ogMs{9,0,t) is superadditive. Therefore, the ratio \g{t) 
decreases monotonically from g{l) = logMc{9) to limt^oo = logm(6*), and Eq. ( |3T] ) 
follows. 

4) We start from [|71 Eq. (7.15)], which states that 

((l-p) , p)(L(e))'QV 

Let m+{9) and m-{9) be the larger and smaller eigenvalues of L(6'), respectively. Inserting 
the spectral decomposition of L{9) with respect to its eigenvalues m+(6*) and m_{9) into 
the expression for Ms{9,0,t), we obtain a constant K{9) such that 

Ms{9,0,t) = K{9){m+{9)y + {l-K{9)){m_{9)Y (Vf > 0, 0 G M). (34) 


Using that Ms{9, 0,1) 


Mc{9), we see that 

My9)-m-{9) 


K{9) = 


m+{9) — m-{9) 
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Since m+{9) > Mc{0) > m_{9), it follows that 0 < K{9) < 1, and we ean drop the 
seeond summand in Eq. ([M|) to obtain Eq. (|3^. 


Remark. The representation for Ms (9) in Eq. (|3^ implies in partieular that 

75(0) = lim ^logM5(0,O,f) = ^logm+( 0 ), 

t^oo Ut U 

in agreement with Eq. ( [27] ). 

The validity of the bounds on the moment-generating funetion and the effeetive capaeity 
extends to a broader elass of time-homogeneous two-state Markov-modulated proeesses, as 
deseribed in (71 Example 7.2.7]. The Markov-modulated proeess is given by a two-state Markov 
ehain in the steady-state and two independent sequenees of i.i.d. random variables, and c\. 
The serviee rate in the fc-th time slot is given by 

Cfc = (1 - Xk)cl + Xkcl , (35) 


where Xk is the state variable of the system. The serviee proeess S{s,t) is defined by S{s,t) = 

E t—1 

k=s ’ 


Lemma 7: Eet S be the two-state Markov-modulated process deseribed above. Assume that the 
random variables c° and c\ have moment-generating funetions and 

If poi +Pio < then the eonelusions of Eemmaj^hold, with Mc{9) = (1 — p)Mco(0) +pMci{9), 
and ... 


L{9) = 


Poo Poi 
PlO Pll 


M,o(0) 0 

0 M,i(0) 


(36) 


Proof: Eix 0 G M, and assume without loss of generality that Mco(0) < Mci{9). We follow 
the proof of Eemma 


Eor the first elaim, we eondition the left hand side of Eq. p9| ) on the state of the system at 
the times ri,..., r„. In a single time slot, we obtain from Eq. (35]) that 


E [e^^'^\Xk] = M,o(0) + (Mei(0) - M,o{9))Xk . 


Sinee the random variables c°, cl and Xk are all independent, it follows that 


E 




Y Y 

y\.Tu . . . -) yXj- 


1[M,o{9) + (M,i(0) -M,o(0))X., 

i=l 


= (^co(0))” '‘"(M,i(0)-M,o(0))''"n^.- 

. 7 c{Ti,...,r„} ieJ 
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Taking expectations, we obtain 


E 


= Y. (M,o(0))" '^'(M,i(0)-M,o(0))'^'e 







Since all coefficients are non-negative, Eq. p9| ) now follows from Lemma 
A similar application of Lemma yields Eq. ( [30l ), which directly implies Eq. ( [3T] ). Einally, 


Eq. follows as before from E Eq. (7.15)]. 


B. Bounds for an MMOO server with feedback 

The results in the previous subsection imply bounds on the equivalent service of a feedback 
system in Eig. containing an MMOO server. We assume d > 0 and te > 0 as parameters of 
the system. Our first result is analogous to Theorem 


Theorem 3: Let S{s,t) be an MMOO service process with feedback. If the transition proba¬ 
bilities satisfy poi +Pio < then, for every 6^ > 0, 

< {{m+i-ey + de-^^)^'^^ , (37) 


where m+{—6) is the larger eigenvalue of the matrix L{—6) defined by Eq. ( [3^ . 

Proof: By Lemma (Parts 1 and 3), we have for each choice of r,,,..., in C„(s, f) the 
bound 


E 


Q-SYJi^iiS{n-i,Ti-d)+s(Tn,t)) 


< Ms{—9, 0,t — s — nd) < {m{—9)y ® . 


Inserting this estimate into the proof of Theorem gives the result. ■ 

The theorem implies the following lower bound on the effective capacity of the service process 
with feedback. 

Corollary 3: Under the assumptions of Theorem the effective capacity 7 „in of the MMOO 
process S{s,t) with feedback satisfies for 6* > 0 

7win(-0) > 7s(-0) - ^ log (1 + , (38) 

where 7 s(—6') is given by Eq. ( |27] ). 

Both results extend to more general two-state Markov-modulated service processes. 

Corollary 4: Let S{s,t) the two-state Markov-modulated process at the end of Subsection IV- 
A. If poi -|- pio < 1, then the moment-generating function of the service process with feedback 
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satisfies Eq. ( [TT] ), where where m+(—6*) is the larger eigenvalue of the matrix L{—6) defined 
by Eq. ( [^ . Its effeetive eapaeity 7 wm(—6^) satisfies Eq. ( [38| ), where 'ysi—0) = —| logm+(—6') 
is the effeetive eapaeity of the serviee proeess. 

Remark: The formula ^{—6) = —|logm+(—6*) agrees with [|71 Eq. (7.17)]. 

Proof: In the proof Theorem replaee Eemma by Eemma |7j ■ 

An important applieation of Corollary is a leftover serviee model with a VBR server (as 


analyzed in See. IV), where eross-traffie arrivals are governed by a two-state Markov-modulated 
proeess. Eet Ck denote the total serviee available to all traffie flows in time slot k, and let al 
denote the eross-traffie arrivals in that time slot. Then the leftover serviee available to the through 
flow ean be bounded from below by 




t-i 


t-i 


^k) 5 


k=s 


- al) 

,k=s 

where we use [x]’*' = max{x, 0}. Corollary applies sinee the sum on the right-hand side is a 
two-state Markov-modulated proeess. It provides a non-trivial lower bound on the serviee of the 
feedbaek system, provided that the stability eondition E[al] < E[ck] is met. 


C. Quality of the MMOO bounds 


As with the VBR server, we ean use the bounds of Theorem to test the aeeuraey of the 
results on the MMOO proeess. Inserting the parameters of the MMOO proeess in Theorem 
S'win is bounded by 


t-i . 

^minjp, < S'„in(s,f) < mini 5 '(s,f), 

k=s ^ 


t — S 


W 


(39) 


The lower bound on the left eorresponds to a serviee proeess of a MMOO proeess with ON 
rate P' = minjP, The lower bound is exaet for when d = 1. Its moment-generating 
funetion ean be bounded with the help of Eemma (Part 3). Eor the upper, we ean only exploit 
Swinisp) < \^]w, sinee we do not have a simple upper bound on the e-quantiles of the 
MMOO proeess. 

We also apply Eemma (Part 4) to the upper bound. Using these expressions in Eq. Q yields 
the following bounds for the effeetive eapaeity: 

Corollary 5: Under the assumptions of Theorems 7win(—6^) is bounded for 9 > Oby 


7win(-^) < 7win(-^) < min|75(-0), , (0 > 0). (40) 
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Here, 74in(~^) the effective capacity of an MMOO process with peak rate P' = min{P, ^}. 
Note that the average rate of this process, given by 7win(0)> hes strictly below the upper bound. 


min{7s'(0, In the limit 6* —?• oo, the three bounds in Eqs. ([3^ and ( |4^ become sharp. 


D. Numerical evaluation of MMOO bounds 

We present numerical examples for an MMOO server with window flow control, proceeding 
in a similar fashion as in the evaluation of the VBR server. The parameters of the MMOO 
process are selected as 

Poo = 0.2, pii = 0.9, P = 1.125 Mb. 


With time slots of length 1 ms, the server has an average rate of 1 Gbps, which is the same rate 
as that of the VBR service evaluated in Subsec. llV-Ci 
In Fig. we present the statistical service curve as a function of time where we 

fix w/d = 100 Mbps in Fig. 8(a) and w/d = 500 Mbps in Fig. 8(b) We set e = 10“®. The 
statistical service curves, plotted as solid lines, have been constructed with the bound on 
from Theorem We include for comparison, upper and lower bounds obtained from Theorem 


via Fq. (j39|) represented by dash-dotted lines. Note that the lower bounds become positive only 
for f > 10 ms, and that this is well matched by the (solid line) statistical service curves. The 
initial latency is a property of the MMOO service process, which may reside for extended time 
periods in the OFF state. As seen in the VBR service in Fig. the statistical service curves 
increase when increasing w and d proportionally. Recall that the lower bound is also a statistical 
service curve when d = 1 ms. As another comparison with the VBR service, we note that the 
upper and lower bounds in Fig. are separated by a wider margin than in Fig. This is due to 
the simpler upper bound, since we have not derived an upper bound for the e-quantiles of the 
MMOO process. 

Fig. 1^ shows the bounds on the effective capacity 7win(—as a function of 6, where we use 
w/d = 100 Mbps in Fig. 9(a) and w/d = 500 Mbps in Fig. 9(b) The bounds from Corollary 
for different values of w and d, are shown as solid lines. The upper and lower bounds obtained 
from Corollary are depicted as dash-dotted lines. When the bounds from Corollary fall below 
those of Corollary the better bound should be used. Except for small values of 9, the lower 
bounds of Corollary are close to the upper bound from Corollary indicating that the effective 
capacity can be accurately computed. 
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(a) w/d = 100 Mbps. 



(b) w/d = 500 Mbps. 

Fig. 8: Statistical service curves 5^jn(0,t) for time-eorrelated (Markov-modulated On-Off) 
service (e = 10“®). 


VI. Conclusions 

We have approaehed a well-known open problem in the stoehastie network ealculus, i.e., an 
extension of the analysis of feedback systems. Our analysis addressed a window flow eontrol 
system with stochastic service, where we considered a serviee with and without time correlations. 
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(a) w/d= 100 Mbps. 



Fig. 9: Effective capacity 7win(—6^) for time-correlated (Markov-modulated On-Off) service. 


We analyzed and then addressed the difficulty of accounting for the time correlations introduced 
by feedback mechanisms. Our analysis revealed major differences between deterministic and 
random window flow control systems, in particular, an additional dependency between the 
characteristic time scales of the feedback delay and the service process. We provided lower 
as well as upper bounds on the available service of the feedback system, which enabled us 
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to discuss the accuracy of our results, and discovered special cases where exact expressions 
for the service can be obtained. The results in this paper can be extended in many directions. 
Obvious generalizations are to consider random feedback delays and time-variable window sizes. 
We chose a window flow control network since a corresponding deterministic network calculus 
analysis exists in the literature. There are numerous other feedback systems that await an analysis 
of their backlog and/or delay properties. Applying our analysis to the detailed dioid algebraic 
models of TCP feedback in [|3| is a logical first step. 
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